import numpy as np
import scipy.linalg as la

n = 80
pmax = 20
w = np.sqrt(list(range(pmax)))

y = np.matrix([
      -3.8709,
      -12.5944,
      32.9643,
      -8.1149,
      25.0225,
      -29.2355,
      -6.3614,
      1.6790,
      -20.9669,
      24.2436,
      16.6332,
      2.3399,
      7.1360,
      -2.5771,
      -14.9766,
      -15.8193,
      -9.4177,
      9.4013,
      11.8224,
      -33.4253,
      -3.0461,
      7.2369,
      0.4587,
      37.6339,
      -22.3718,
      27.0435,
      -11.8499,
      0.9562,
      -0.7392,
      -11.6313,
      37.8674,
      4.9134,
      -26.1744,
      -27.3233,
      -7.4307,
      -20.7317,
      5.0787,
      20.2819,
      17.0273,
      -10.6180,
      7.2649,
      9.5919,
      -10.5182,
      -3.7782,
      26.4167,
      -18.3995,
      -5.4334,
      -14.2456,
      11.1707,
      13.6572,
      -3.5610,
      25.8739,
      -17.8474,
      -1.4456,
      -33.4935,
      -8.9574,
      18.9878,
      -0.5333,
      -15.5718,
      17.9230,
      2.0426,
      -9.6112,
      6.1276,
      -1.9454,
      16.7150,
      5.6382,
      8.7873,
      4.3659,
      -22.2739,
      29.2610,
      -11.5744,
      -39.7351,
      10.8082,
      -23.9951,
      -2.0934,
      6.9986,
      5.2206,
      -0.3450,
      12.0607,
      26.6712 ]).T
